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1  Introduction 


This  project  addressed  the  interaction  between  instability  waves  and  Klebanoff 
modes  for  their  relevance  to  compressors.  It  provided  peripheral  support  to 
computer  simulations  of  transition  on  compressor  blades  (see  §5:  Papers  com¬ 
ing  from  this  work).  In  light  of  the  re-orientation  of  AFOSR  university  pro¬ 
grams,  I  should  note  that  the  proposal  was  funded  under  Dr.  Jeffries’  program 
on  turbomachinery  and  rotating  flow. 

It  is  generally  believed  that  in  turbomachines  transition  is  via  bypass  mecha¬ 
nisms.  There  was  a  hint  in  the  literature  that  on  compressors  instability  waves 
might  play  a  role.  But  that  came  from  a  good  deal  of  filtering  and  wavelet 
transforming  of  a  few  transducer  signals  on  a  compressor  rig.  So  we  proposed 
to  explore  the  topic  at  a  basic  level,  under  this  AFOSR  grant.  We  were  quite 
surprised  to  find  that  the  interaction  between  Klebanoff  modes  (the  precur¬ 
sors  to  bypass  transition)  and  Tollmein-Schichting  waves  (which  are  what  is 
bypassed)  could  cause  abrupt  transition:  discussed  the  following  sections  of 
this  report.  About  the  time  of  this  grant  we  were  initiating  DNS  of  transition 
on  a  compressor  geometry  (most  of  the  funding  was  via  collaboration  with 
colleagues  at  Univ.  of  Karlsruhe  in  Germany.  Those  simulations  were  carried 
over  a  couple  of  years  and  a  last  paper  is  only  now  under  review) .  They  confirm 
that  the  Klebanoff  mode-TS  interaction  does  indeed  take  place  on  compressor 
blades.  This  is  contrary  to  the  situation  on  turbine  blades,  and  probably  is 
influenced  by  the  adverse  pressure  across  a  compressor. 

The  following  summarizes  some  of  the  basic  research  done  under  this  program. 
It  is  motivated,  after  the  fact  by  the  image  in  figure  1  from  our  DNS  of 
transition  on  a  compressor  blade,  beneath  free-stream  turbulence.  That  figure 
is  showing  a  transition  mechanism,  in  context.  We  started  by  studying  the 
potential  for  such  a  process  in  the  abstract,  under  AFOSR  support. 


Streamwise  length 

Fig.  1.  Plane  and  side  views  of  the  tangential  perturbation  field  on  DNS  of  a 
compressor  blade.  The  dashed  line  marks  the  location  of  the  side  view  (bottom 
pane).  The  side-view  plane  bisects  a  disturbance  wave  which  leads  to  breakdown  to 
turbulence. 

2  Simulations  of  the  basic  study 


The  numerical  method  for  this  study  is  a  finite  volume,  staggered  mesh,  frac¬ 
tional  step  algorithm  for  the  time-dependent,  three-dimensional,  incompress¬ 
ible  Navier-Stokes  equations  in  generalized  coordinates.  The  spanwise  direc¬ 
tion  is  assumed  to  be  periodic  and  is  treated  by  spectral  methods  to  reduce 
computational  cost.  The  convective  terms  and  the  off-diagonal  diffusion  terms 
are  advanced  explicitly  by  the  Adams-Bashforth  scheme,  while  the  remaining 
diffusion  terms  are  treated  implicitly  by  Crank-Nicolson.  The  Poisson  equa¬ 
tion  for  pressure  is  solved  by  line  relaxation  with  multi-grid.  The  code  is  MPI 
parallel  and  simulations  were  performed  on  32  processors  of  an  AMD  Opteron 
280  cluster. 

The  inflow  is  specified  as  a  superposition  of  a  3-D,  spatially  evolving  continu¬ 
ous  mode  (uCOn)  and  a  2-D  spatially  unstable  Tollmien-Schlichting  wave  (uts) 
on  a  Blasius  mean  flow  ( ): 


uq  =  Ub  +  Real 
vo  =  Vb  +  Real 


Wq  =  Real 


Ar 


[■ Aconucon(y)e~WQ3t  cos (kzz)  +  ATSuTs(y)e 
\AconVcon{y)e~lUost)  cos (kzz)  +  ATSvTs{y)e 
wC(m{y)e~lUost)  sin (kzz) 
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The  Tollmien-Schlichting  and  continuous  modes  are  obtained  by  solving  the 
Orr-Sommerfeld  and  Klebanoff  equations  with  well  established  methods:  a 
Chebyshev  collocation  scheme  is  used  to  find  the  discrete  modes  and  an  im¬ 
plicit,  matrix  method  is  used  for  the  continuous  modes.  The  inflow  T-S  wave 
has  a  non-dimensional  frequency 

LUIS  R 

F  S  Jp10  =  124- 
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At  the  inlet  Reynolds  number,  Rsgg  =  2, 000,  it  is  unstable  and  has  a  complex 
wavenumber 

ch599  =  0.6643  -  3.355  x  10“3i.  (2.2) 

The  inflow  continuous  O-S  modes  have  non-dimensional  frequency,  wall-normal, 
and  spanwise  wavenumbers  of 

F  =  33,  kySgg  =  7r/2,  kz5gg  =  n2iT  /  Lz  =  nkQz  (2-3) 
The  spanwise  ( z )  wavenumber  is  expressed  as  a  multiple  of  a  wave  number 

k°  =  y:  (2-4) 

having  period  equal  to  the  domain  width.  The  domain  width  is  eight  times 
the  initial  boundary  thickness,  Lz  =  8^99.  n  is  the  number  of  waves  spanning 
the  domain.  Hence,  the  physical  scale  of  an  n-wave  inlet  disturbance  is  8Sgg/n. 
Results  will  be  presented  for  n  =  2  and  n  =  5  —  which  will  be  called  mode  2, 
or  long  wavelength,  and  mode  5,  or  short  wavelength.  It  should  be  emphasized 
that  tests  with  a  wider  domain,  and  with  these  same  wavelength  disturbances, 
demonstrated  domain  independence.  The  relevant  quantity  is  the  wavelength 
86gg/n,  not  the  number  of  waves  in  the  computational  domain. 
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Fig.  2.  Contours  of  fluctuating  velocity  of  each  mode  alone:  top  pane,  T-S  wave; 
middle  pane,  continuous  mode.  Both  modes  are  present  in  the  lowest  pane. 
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3  Results 


Either  mode,  of  itself,  is  unable  to  provoke  transition.  The  two  top  panes  of 
figure  2  contain  contour  plots  of  the  T-S  wave  and  the  continuous  mode  2 
alone.  The  T-S  wave  grows  exponentially,  but  remains  two-dimensional.  The 
continuous  mode  is  three-dimensional,  but  stable.  It  induces  a  disturbance  in 
the  boundary  layer,  that  evolves  into  a  jet-like  form:  the  velocity  contours  are 
elongated  in  the  x  direction  and  the  dominant  velocity  component  is  u.  Be¬ 
neath  free-stream  turbulence,  these  jets  are  often  termed  streaks,  or  Klebanoff 
‘modes’,  or  better,  Klebanoff  ‘distortions’.  The  growth  of  streamwise  elongated 
disturbances  can  be  understood  via  linear  theory.  However,  quadratic  nonlin¬ 
earity  plays  a  role  in  the  middle  case  of  figure  2.  The  spanwise  spacing  of  the 
contours  is  halved  toward  the  end  of  the  domain.  This  is  caused  by  quadratic 
non-linearity  acting  on  the  initial  spanwise  periodicity.  However,  it  is  seen  in 
the  bottom  pane  of  this  figure  that  transition  occurs  before  the  wavelength 
has  been  halved. 

With  both  modes  present  transition  usually  occurred  within  the  computational 
domain.  That  is  illustrated  in  the  lowest  plot  of  figure  2  and  by  figure  3.  The 
later  are  skin  friction  curves  comparing  simulations  with  the  unstable  T-S  wave 
alone,  and  with  both  modes  present  at  the  inlet.  Lines  presenting  Cf  versus 
Rx  in  laminar  and  in  turbulent  flow  are  included  for  reference.  The  average 
of  the  T-S  wave  alone  case  falls  on  top  of  the  laminar  line.  With  both  modes 
present,  the  flow  transitions,  overshooting  the  turbulent  line.  That  overshoot 
is  typical  of  transition,  although  it  can  be  exacerbated  if  the  high  demand  for 
streamwise  grid  resolution  once  the  boundary  layer  becomes  turbulent  is  not 
fully  met. 


Fig.  3.  The  skin  friction  of  modes  2  and  5.  The  discrete  and  continuous  modes  both 
have  a  1%  amplitude  at  the  inlet. 

The  pattern  of  the  perturbation  when  transition  occurs  is  illustrated  at  the 
bottom  of  figure  2.  Transition  is  preceded  by  the  appearance  of  A  shaped 
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velocity  contours.  Although  this  is  reminiscent  of  secondary  instability  of  Toll- 
mien- Schlichting  waves,  the  lateral  spacing  between  A’s  is  very  much  narrower 
and  seems  to  be  controlled  by  the  spanwise  mode  number,  n,  of  the  continuous 
mode.  Two  broad  classes  of  behavior  were  seen,  as  epitomized  by  modes  2  and 
5  (c.f.,  equation  2.3).  Those  two  cases  are  discussed  in  the  next  sections. 

3.1  Long  spanwise  wavelength 


A  simulation  with  only  the  longer  wavelength  mode  prescribed  at  the  inlet 
is  summarized  in  figure.  The  horizontal  coordinate  is  measured  relative  to 
the  inlet  at  Xq  ==  80.  Figure  4  looks  down  on  a  horizontal  plane.  The  inlet 
Orr-Sommerfeld  disturbance  has  the  spanwise  wavelength  and  streamwise  pe¬ 
riodicity  of  the  continuous  mode.  Its  non-linear  evolution  is  illustrated  by  u 
contours. 
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Fig.  4.  Contours  of  streamwise  fluctuating  velocity  in  a  horizontal  cross  sectional 
plane:  mode  2. 

The  vertical  sections  in  figure  5  show  the  Klebanoff  response  more  clearly.  The 
Orr-Sommerfeld  mode  is  seen  in  the  free-stream.  Klebanoff  modes  are  created 
inside  the  boundary  layer.  They  distort  the  disturbance.  Initially,  the  phase 
varies  with  height  inside  the  boundary  layer.  The  disturbance  subsequently 
intensifies  —  each  plane  is  scaled  between  its  local  maximum  and  minimum  so 
that  contours  in  the  free  stream  become  gray.  Nonlinearity  becomes  apparent 
by  x  —  xq  =  80,  as  the  spanwise  wavelength  begins  to  halve  next  to  the 
wall.  The  cross-sectional  planes  show  how  the  wavelength  doubles  from  the 
wall,  outward  to  the  free-stream:  near  the  top  of  the  boundary  layer  the  dark 
contours  are  seen  to  merge,  forming  an  inverted  ‘v\ 

By  x  —  Xo  =  170  the  boundary  layer  disturbance  has  become  quite  distinct 
from  that  in  the  free-stream.  The  lower  image  in  figure  5  shows  the  variation 
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Fig.  5.  Contours  of  streamwise  fluctuating  velocity  in  vertical  cross  sectional  planes. 
The  plane  at  x  —  xq  =  180  is  shown  at  the  bottom. 

of  spanwise  wavelength  with  height  quite  clearly. 

Returning  to  figure  4,  the  wavenumber  doubling  is  seen  in  the  z  spacing  of  the 
dark  contours.  The  surface  elevations  in  the  lower  part  of  the  figure  show  how 
undulations  in  the  x  direction  slowly  decay,  leaving  a  streamwise  elongated 
disturbance.  When  the  T-S  wave  is  added  to  mode  2  at  the  inlet,  transition 
occurs  before  the  spanwise  wavelength  is  halved,  as  will  be  seen  shortly.  The 
evolution  of  mode  2  far  downstream  is  superseded  by  transition. 

With  both  the  continuous  and  the  discrete  mode  present,  transition  occurs. 
Upstream  of  full  transition  to  turbulence,  velocity  contours  develop  a  A  shaped 
pattern.  The  A’s  have  a  spanwise  wavelength  equal  to  that  of  the  continuous 
mode. 

Figure  6  is  an  instantaneous  view  of  the  perturbation  field.  It  captures  the 
early  stages  of  development  of  a  pair  of  A-structures  at  x  ~  30,  and  a  mature 
pair  further  downstream,  at  x  ~  38.  The  A’s  are  shaded  by  the  mean  velocity. 
Thereby,  lighter  gray  indicates  greater  distance  from  the  wall.  The  shaded 
surface  is  defined  by  Q  =  —0.01,  where  Q  =  diUjdjUi  is  the  difference  between 
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Fig.  6.  Instantaneous  disturbance  showing  the  early  stage  of  a  pair  of  A  structures, 
second  pair  downstream  and  breakdown  to  turbulence.  Constant  Q  surface,  shaded 
by  mean  velocity  over  velocity  vectors  in  a  plane. 

rate  of  strain  squared  and  rate  of  rotation  squared.  A  region  with  Q  <  0  is 
considered  to  identify  a  vortex.  The  constant  Q  surfaces  show  vortices  lifting 
up,  then  breaking  down  into  turbulence.  Velocity  vectors  show  both  the  jets 
caused  by  the  continuous  modes,  and  the  flow  induced  by  the  vortices.  These 
aspects  will  be  pursued  in  subsequent  figures. 

90_ 95_ 1^0 


Fig.  7.  Time  sequence  showing  the  inception  of  the  A-structure,  its  growth,  and 
finally  breakdown.  Smooth  surfaces  are  the  set  Q  =  —  0.01. 


Figure  7  is  a  time  sequence  of  the  inception,  growth,  and  breakdown  of  A- 
structures.  Note  that  the  viewing  window  translates  downstream  from  frame 
to  frame.  Vectors  show  the  in-plane  perturbation  velocity  field.  Light  arrows 
represent  forward  velocities  and  darker  shades  indicate  velocities  that  are  neg¬ 
ative  relative  to  the  mean  flow.  Constant  Q  surfaces  cross  the  plane  of  velocity 
vectors.  As  previously,  they  are  shaded  by  the  mean  velocity,  so  dark  regions 
are  near  the  wall. 
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Proceeding  from  the  upper  left,  downward  in  the  first  snapshot  a  vortical 
region  is  observed  near  x  =  15.  In  the  second  snapshot,  this  region  becomes  the 
inception  point  for  a  small  vortex  leg,  which  commences  on  a  positive  velocity 
perturbation  (light  vectors).  A  second  leg  emerges  in  the  third  snapshot,  atop 
a  region  of  negative  velocity  vectors.  The  gray  scale  shading  shows  that  the 
vortices  are  lifting  away  from  the  wall,  angled  up  from  left  to  right. 


The  lightest  shade  of  gray  is  near  the  outer  edge  of  the  boundary  layer.  In 
subsequent  snapshots,  the  A  vortices  extend  increasingly  upward,  reaching  the 
top  of  the  boundary  layer  as  they  convect  downstream. 


Fig.  8.  The  breakdown  of  the  A  structures.  Time  increases  from  upper  left  to  lower 
right. 

A  clear  asymmetry  is  captured  in  the  time-sequence  and  is  unique  to  this  tran¬ 
sition  scenario  (unlike  the  secondary  instability  of  T-S  waves  in  the  absence 
of  streaks).  Despite  emerging  first,  the  lower  leg,  along  the  forward  streak,  is 
appreciably  shorter  than  the  upper  leg.  In  the  field  of  velocity  vectors,  vortical 
motion  associated  with  the  A’s  is  only  clear  around  the  lower  leg.  The  upper 
leg  is  significantly  elongated  due  to  the  backward  streak,  which  dominates  the 
perturbation  field  in  the  plot  of  velocity  vectors.  The  A’s  continue  to  stretch 
and  lift  away  from  the  wall,  and  finally  breakdown  to  turbulence  at  the  lower 
right. 

Figure  8  is  a  time  sequence  showing  contours  of  w-perturbations  in  a  cross¬ 
stream  plane.  The  time  instance  in  the  top  left  pane  is  an  early  stage  in 
the  evolution  of  the  A-structures.  The  cross-sectional  plane  is  located  down¬ 
stream  of  the  emerging  vortices.  The  u  velocity  contours  demonstrate  that 
the  boundary  layer  disturbance  is,  at  this  instant,  still  at  the  fundamental 
spanwise  wavenumber.  The  incident  continuous  mode  is  apparent  in  the  free- 
stream.  Hence  the  wavenumber  doubling,  seen  in  figure  2,  would  occur  further 
downstream,  and  is  superseded  by  the  modal  interactions  and  transition  to 
turbulence. 

At  all  subsequent  times,  the  cross-sectional  plane  is  traversing  the  A’s.  A  region 
of  negative  w-perturbation  (dark  contours)  is  observed  between  the  two  legs 
of  the  structure,  and  is  contained  within  an  originally  forward  (light  contours) 
streak.  This  negative  u-perturbation  is  induced  by  the  vortical  structure  itself. 


Circulation  round  the  A  shape  causes  upwelling  in  between  the  legs.  This 
carries  low- momentum  fluid  from  the  near  the  wall.  Perhaps,  it  should  be 
noted  that,  beyond  the  cross-sectional  plane,  the  spanwise  wavenumber  of  the 
continuous  mode  can  be  observed,  even  after  breakdown  to  turbulence.  Bands 
of  turbulence  emanate  from  the  loop  at  the  top  of  the  A  vortices.  When  the 
boundary  layer  becomes  fully  turbulent,  the  continuous  mode  wavelength  is 
no  longer  evident  (figure  2). 
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